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Recently a procedure by generalized density matrix (GDM) is proposed [lj for calculating a 
collective/bosonic Hamiltonian microscopically from the shell- model Hamiltonian. In this work we 
examine the validity of the method by comparing the GDM results with that of the exact shell- 
model diagonalization in a number of models. It is shown that the GDM method reproduces the 
low-lying collective states quite well, both for energies and transition rates, across the whole region 
going from vibrational to 7-unstable and deformed nuclei. 
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It is a long-standing problem in nuclear physics to 
understand how macroscopic collective motion arises 
from microscopic single-particle motion. The shell-model 
(configuration interaction) successfully reproduces var- 
ious collective behaviors by diagonalizing the nucleon 
Hamiltonian in a huge Slater-determinant basis. How- 
ever, the dimension of the basis makes it impractical ex- 
cept for the cases with only a few valence nucleons. On 
the other hand, phenomenological bosonic approaches are 
often successful in fitting the experimental data (first of 
all the geometric Bohr Hamiltonian 2, 3] and the inter- 
acting boson model [4]). This shows that, out of the 
huge Slater-determinant space, there exists a few de- 
grees of freedom with a bosonic nature, which are usually 
enough in describing the collective states. Serious efforts 
were devoted to deriving those parameters of the bosonic 
Hamiltonian from the underlying shell-model Hamilto- 
nian. However, the complete theory is still missing. 

Recently we proposed |l[ a procedure based on the 
generalized density matrix (GDM) that was originally 
formulated in Refs. This procedure is rather sim- 

ple, clean, and consistent. In compact form, there are 
only two equations, (14) and (23) in Ref. [l(. Results 
from the lowest orders give the well-known Hatree-Fock 
(HF) equations and random phase approximation (RPA). 
Higher orders fix the anharmonic terms in the collec- 
tive/bosonic Hamiltonian. The aim of this work is to 
demonstrate the validity of the GDM method, by com- 
paring its results with that of the exact shell-model di- 
agonalization. 

In this work for simplicity we restrict ourselves to sys- 
tems without rotational symmetry. The GDM formula- 
tion with angular-momentum vector coupling has been 
considered in Ref [9(. The single particle (s.p.) space in 
this work is drawn schematically in Fig. [T] There are two 
degenerate s.p. levels with energies e = ±1/2. The Fermi 
surface is in between, thus the lower levels arc completely 



filled and upper levels are empty. Each s.p. level has a 
quantum number m that is a half integer. Degenerate 
time-reversal pair has m with different sign, toj = —mi. 

For fermions, |1) = — and we choose the phases such 
that 

\fh) = I — m) , I — m) = — |m) (to > 0). 
We assume a two-body Hamiltonian, 

H = ^ fi 2 a\a 2 + 7 ^ V 12 34N[a\ala 3 a4\, (1) 

12 1234 

where fi 2 = S12S1, e\ are the HF s.p. energies shown 
in Fig. [TJ The density matrix p\ 2 = S± 2 ni, where the 
occupation number n\ = 1(0) for the lower(upper) s.p. 
levels. ./V[a{a2a3a4] is the normal-ordering form of oper- 
ators. The residual interaction is of Q ■ Q type, 

Vl234 = -K(gi4<723 - 913924), 

where the moment operator Q = ^2 12 qi2a\ a 2 is Hermi- 
tian and time-even. For simplicity we assume q is real, 
thus 

912 = 921 = 921 = 9l2- 

Operator q has certain selection rule with respect to 
quantum number m, which will be specified later. We 
further set diagonal matrix elements of q to zero, 911 = 0; 
hence the moment of the mean field vanishes, Q(°°) = 
Tr[9p] =0. 

Following the procedure in Ref. |l| , we are able to map 
the fermionic Hamiltonian ([T]) onto a bosonic Hamilto- 
nian 

Hb = J± + *L + A(30)^! + A (12) {Jhfl + A (40)4 

2 2 3! 4 4! 

+A (22)KV} +A (04)^ +A (50)^! + _ (2) 

8 4! 5! 



'Electronic address: jial@nscl.msu.edu 



where the collective coordinate a and momentum it sat- 
isfy [a, 7r] = i. Meanwhile, the generalized density matrix 
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r (mn) j n fog expansion 
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D — t , (10) , (01) , (20) a ,„■. 

i?i2 = a 2 ai = pi2 + r{ 2 'a + r\ 2 >n + r\ 2 — + . . . (3) 

is expressed in terms of K^ mn \ The bosonic Hamiltonian 
Hf, ([2} should reproduce the low-lying collective spectrum 
of the original nucleon Hamiltonian H (fl]) . Substituting 
the solution ([3]) into Q = J2i2 qi2a\a,2, we get the boson 
image of the latter, 

Q b = Q™a + + Q(° 2 >y + + ..., (4) 

where Q( m ") = Tr[qr( mn )], and time-odd terms vanish 
automatically. The transition rates calculated from Qb 
between eigenstates of Hb should reproduce that of Q 
between eigenstates of H. 

As stated in Ref. [1], the GDM method fixes H b 
completely. In each even order (quadratic, quartic 
...) in Hi,, the GDM method gives one constraint on 
A(" m )'s. The number of constraints is the same as 
that of independent parameters in Hb, removing in 
Eq. © superficial degrees of freedom due to canonical 
transformations of a and it conserving [a, 7r] = i. 

In the following we consider four models with different 
structures (different configurations of s.p. levels and dif- 
ferent selection rules of q). We start with the simplest 
case of the first model. Both the upper and lower group 
have 12 degenerate s.p. levels with quantum numbers 
m = ±5,±|,...,±Tp Operator q has the selection rule 
Am = 0, that is, q±2 vanishes unless mi — m,2- The 
non-vanishing q\ 2 (mi = m 2 ) are set to be 1. 

In this model we find by numerical computation a 
"symmetry" . That is, there are only three non- vanishing 
terms in Hamiltonian ([2]): uj 2 , A^ 40 \ and A^ 22 ) (besides 
7r 2 /2). We suspect there was some kind of "quasi angular- 
momentum symmetry", similar to that in the Lipkin 
model, where the only three non- vanishing terms are ui 2 , 
A( 4 °), and A(° 4 ) (see Ref. fl). 

The results are shown in Fig. [5] We see that the 
GDM calculation reproduces the exact results of the shell 
model quite well, both for energies and transition rates. 
In the shell model we calculate the lowest several states 
by Lanczos method. The dashed line in the upper panel 
is the beginning of the s.p. continuum, only those col- 
lective states below the continuum can be seen. In the 
GDM calculation the resultant bosonic Hamiltonian is 
diagonalized in a finite "physical" bosonic space, {|0 < 
n < 12)}. \n) is the n phonon state, j47A|n) = n|n), 
A* = (ua + W7r)/v2, uv = — 1. The coefficient u is fixed 
by minimizing A|HF) in its one-particle-one- hole compo- 
nents, where |HF) is the Hartree-Fock ground state that 
is mapped onto the bosonic state \n = 0). The result 
T k (10) | 2 

is u — <F<1 mn , where the summation indices 1 

Z^2<F<1 I' 12 I 

and 2 run over unoccupied and occupied s.p. levels, re- 
spectively ( ".F" means Fermi surface) . In models of this 



work, u is a number close to 1. The shown GDM ener- 
gies and transitions are practically independent of small 
variations of u around the above value. 

As k increases, the system goes from vibrational to 
7-unstable region. In the vibrational region, higher ex- 
cited states are influenced more by the anharmonicitics, 
as expected. At large k the spectrum becomes doubly de- 
generate in a deep double-well potential (large negative 

We would like to mention an important point, that the 
GDM method works better with increasing ft (collectiv- 
ity). Another calculation has been done (but not shown 
here) with 8 particles in 16 s.p. levels. The GDM results 
of the current calculation (12 particles in 24 s.p. levels) 
have very clear improvement over that of the former. In 
other words, the error in Fig. [5] may be of order 1/51. 
The largest part of this error may come from the RPA 
frequency ui 2 . At the current stage, the GDM method 
calculates all A^ mn > in their leading order of but 
not the next. We suspect that the correct uj 2 = A' 20 ) is 
smaller by a term than the one determined here by 
the RPA equation. This would shift all the GDM curves 
to the left (smaller n) a little bit, which would decrease 
greatly the systematic error (see Fig. [2]). This system- 
atic error own to uj 2 seems to be present in all models 
in this work (see Figs. [3][5]). Also, it is confirmed in 
the Lipkin model where everything is known analytically 
(see Ref. [l|). Hence an apparent improvement would be 
calculating uj 2 in its next-to-leading order of I/O. 

Next we consider the second model, which has the same 
s.p. configuration but a different q that now has the 
selection rule Am = 0,±1. Non-vanishing qi2 is still set 
to be 1. In this model we did not find a "symmetry" as in 
the previous case, so the problem exists of what should be 
the "best" mapping. In the following we did three sets 
of GDM calculations. The first calculation keeps only 
A( 40 ) (besides w 2 a 2 /2 and 7r 2 /2) in Hb, which is fixed 
by the constraint from 4th order in e.o.m. The second 
calculation keeps only A^ 40 ) and A' 60 ', which are fixed by 
the two constraints from upto 6th order in e.o.m. The 
third calculation keeps only A (40) , A^ 22 ), and A (04) , fixed 
by the three constraints from upto 8th order in e.o.m. 

We first notice in Fig. [3] that in this model the s.p. 
continuum goes down as increasing k, as opposed to going 
up in the previous model. This is because now mixing 
of s.p. levels within the upper(lower) group is allowed 
by the selection rule that Am can be ±1. As a result, 
originally degenerate levels from the upper (lower) group 
get a finite spread, which decreases the gap of the s.p. 
continuum. Only the first excited state is below the gap 
and calculated in the shell model. 

In the GDM calculations we see that the simplest one 
degree-of- freedom (A^ 40 )) calculation is reasonably well 
in most cases except at very large k. The other two cal- 
culations (A( 4 °/ 6 °) and A( 40 / 22 / 04 )) give essentially the 
same results (for the quantities shown in Fig. [3]), al- 
though their common parameter, A^ 40 \ differ a lot. This 
insensitivity of GDM results to the degrees of freedom 
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chosen, is important. As we said, two different bosonic 
Hamiltonians could be equivalent if they were related by 
canonical transformations/renormalizations of variables 
a and tt. This insensitivity simply says that the GDM 
formulism knows these renormalizations and does them 
correctly. In the first model we also hnd this insensitiv- 
ity (but not shown). Finally we notice that in regions of 
bj 1 ~ I/O, calculations that go to higher orders in e.o.m. 
may give unphysical results. This is again because the 
e.o.m. are accurate in leading order of l/f2 but not the 
next. The fact that this "divergence" appears slightly 
before the instability point of RPA shown in Fig. [3l indi- 
cates again that the correct uJ 2 may be smaller than the 
one calculated by RPA, as mentioned before. 

At last we consider two models with s.p. configura- 
tion that is asymmetric in upper and lower groups, which 
generates odd anharmonicities that are necessary for de- 
formation. In the third model, the lower group has 10 
s.p. levels with m = ±|, ±4p the upper group has 14 
s.p. levels with m = ±i, ±4f-. In the fourth model, 
the lower group has 12 s.p. levels with m = ±|, 
the upper group has 10 s.p. levels with m = ±i, ...,±|. 
In both models, operator q still has the selection rule 
of Am = 0, ±1, with non- vanishing matrix elements set 
to be 1. The third model has a slightly larger asymme- 
try than that of the fourth model, and their signs of the 
asymmetry are different. 

These two models are more complicated in the sense 
that now there are more active degrees of freedom (odd 
anharmonicities). In the GDM method, we do a possi- 
bly simplest calculation. We keep in Ht only A^ 30 ), A' 12 ', 
and A( 4 °) (besides u 2 a 2 /2 and tt 2 /2). A( 30 ) and A< 12 ) are 
fixed by requiring Q( 20 ) = = in the solution (0}. 

Then A^ 40 ) is fixed by the constraint from 4th order in 
e.o.m. The requirement = = is the same as 

that for previous two models without upper-lower asym- 
metry, by which A^ 30 ) and A^ 12 ) vanish. 

The results are shown in Figs. 0][5] The deformation 
begins around the critical point of RPA when oJ 2 becomes 
negative. In the vibrational region the potential is stiff 
and deformation is not easy. As k increases, the poten- 
tial becomes flat in bottom and finally of a double-well 
shape. Then, even a relatively small odd anharmonicity 
(here mainly A^ 30 )) can tilt the potential and generate 
large deformation. We notice firstly that the GDM 
calculations give the correct sign of deformations: 
positive/negative for the ground/first-excited state of 
the third model, and vice versa for the fourth model. In 
realistic situation A( 30 )((a x a) 2 x a) ~ A^ 30 )/? 3 cos 37 
(a is the quadrupole phonon and /3, 7 are Bohr angles), 
the sign of A^ 30 ) "determines" the intrinsic shape of the 
nucleus (prolate or oblate). This is especially interesting 
in the transitional regions where the rotor formula is 
not applicable. Secondly, the quantitative agreement 
of deformation is also good except at the largest k. 



There the deformation "saturates" towards its maximal 
possible value within the model space, favored by energy. 
Meanwhile in the boson mapping, we are too close to 
the boundary of the finite physical bosonic space, and 
the GDM results become inaccurate. In realistic nuclei 
this "saturation" may not be usual. The number of 
participating/active nucleons is usually around 30 in 
well-deformed nuclei, which is much larger than that 
around 10 in the current models. Finally, we would like 
to point out that the first excited state in our simple 
models is not a "rotational" state, rather it corresponds 
to the next "band head" in realistic rotational nuclei. 
The rotational states that are very low in energy come 
in only in three dimensions. 

In this work we demonstrate the validity of the 
GDM procedure in microscopically calculating the col- 
lective/bosonic Hamiltonian. The lowest several states 
of this bosonic Hamiltonian quite well reproduce the col- 
lective states of the exact shell-model, both for energies 
and transition rates, in a wide range from vibrational, 7- 
unstable, to deformed nuclei. Specifically, we show that 
deformation can be described without introducing a de- 
formed mean field. The traditional procedure of "sym- 
metry breaking and restoration" , first "statically" breaks 
rotational symmetry in the ground state, by representing 
the latter as a Slater determinant of deformed s.p. levels 
(Nilsson levels); then projects afterwards to good angular 
momentum. However, in case of large shape fluctuations 
(flat minimal of energy surface) or shape coexistence (two 
close minimal), it may fail. On the other hand, the GDM 
procedure always conserves the rotational symmetry. De- 
formations are put in "dynamically" at higher orders (for 
example cubic terms) beyond the mean field. Thus it is 
suitable to describe for example shape fluctuation and 
coexistence problems. 

In realistic nuclei, the gap of s.p. continuum is gen- 
erated by the pairing correlations. The GDM formulism 
based on the Hartree-Fock-Bogoliubov transformation is 
straightforward and has been done in Ref. [9]. How- 
ever, another treatment may be possible. Instead of 
introducing Bogoliubov quasi-particles and representing 
the ground state as a Slater determinant of the former, 
the pairing correlations are considered in higher orders 
beyond the mean field, by keeping both the particle- 
hole and particle-particle channels in the factorization 
a\ala,2ai w a\a\ ■ a\a,2 ~ a\a,2 • a\a\ + o\al ■ a2&i. In this 
way the symmetry of exact particle number is always 
conserved. Work along this line is in progress and results 
seem promising. We are also generalizing the GDM code 
by including angular-momentum vector coupling that are 
necessary for realistic calculations. 

This work is supported by the NSF grants PHY- 
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FIG. 1: Single-particle level scheme. 1 is the time-reversal level of 1. Each level has a quantum number 
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FIG. 2: (Color online) Excitation energies E n — Eo and transition matrix elements (m\Q\n) in the first model of this work as 
a function of n. The black solid lines are exact results by shell-model diagonalization. The black dashed line is the beginning 
of the single particle continuum. The red dashed-dotted lines are the RPA results. The blue symbols are GDM results. The 
asterisks are energies; and the squares, circles, up-triangies and down-triangles are matrix elements of Q between different 
states. (m\Q\n) that are not shown vanish in both the shell model and the GDM calculations. 
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FIG. 3: (Color online) Excitation energies E„ — Eo and transition matrix elements (0|Q|1) in the second model of this work as 
a function of k. The black solid lines are exact results by shell-model diagonalization. The black dashed line is the beginning 
of the single particle continuum. The red dashed-dotted lines are the RPA results. The blue symbols (asterisks, up-triangles 
and circles) are results of three different sets of GDM calculations, as labeled in the legend. 
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FIG. 4: (Color online) Excitation energies E n — Eq and transition matrix elements (m\Q\n) in the third model of this work as 
a function of k. The black solid lines are exact results by shell-model diagonalization. The black dashed line is the beginning 
of the single particle continuum. The red dashed-dotted lines are the RPA results. The blue symbols are GDM results. The 
asterisks are energies; and the squares, circles and up-triangles are matrix elements of Q between different states. 
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FIG. 5: (Color online) Excitation energies E n — Eo and transition matrix elements (m|Q|n) in the fourth model of this work as 
a function of k. The black solid lines are exact results by shell-model diagonalization. The black dashed line is the beginning 
of the single particle continuum. The red dashed-dotted lines are the RPA results. The blue symbols are GDM results. The 
asterisks are energies; and the squares, circles and up-triangles are matrix elements of Q between different states. 



